Detect AG at BNL-RHIC via Double Quarkonium Production 
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The double spin asymmetry for exclusive J/^l> pair production in the polarized p — p collisions at 
RHIC is investigated. Our study shows that the asymmetry in this process is measurable at RHIC- 
SPIN experiments in the near future, and, hence it can be used to determine the gluon distribution 
function /S.G(x) in the polarized nucleon. 
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Since the first measurement of the polarized structure 
function by the European Muon Collaboration [1], an 
enormous amount of researches have been done on the 
nucleon spin structure both experimentally and theoret- 
\ ically. The unpolarized deep inelastic scattering (DIS) 
> experiments show that the gluon shares a large portion 
of the parent proton's momentum. However, how the 
gluons share the spin of proton is still an open question. 
' In the conventional DIS process, the gluon contributions 
are the Quantum Chromodynamics(QCD) higher order 
[ effects and we can not avoid some ambiguities coming 
• from e.g. the factorization scheme dependence, when de- 
' termining the polarized gluon distributions. Although 
there are some efforts [2,3] to parameterize them, it seems 
difficult to obtain clear understanding of the role of glu- 
ons. 

Therefore, it is desirable to study other processes than 
DIS in which a direct measurement of the polarized gluon 
. is possible. It is now expected that the polarized proton- 
' proton collisions (RHIC-Spin) at BNL relativistic heavy- 
, ion collider RHIC will provide copious experimental data 
' to unveil the polarized parton distributions. Since the 
, major emphasis and strength of RHIC-Spin is to mea- 
' sure the gluon polarization, it is important and interest- 
ing to investigate various processes which are attainable 
experimentally to this aim. 

Up to now, three processes are considered to be promis- 
ing for measuring the polarized gluons, which are thought 
to be feasible at RHIC technically. Those are 

• High-pT Prompt Photon Production 

• Jet production 

• Heavy Flavor Production 

There are some advantages and disadvantages in each of 
these processes. The prompt photon production with po- 



larized beams at RHIC is a useful process, but we must 
know precisely the quark distributions from the polarized 
DIS experiments. For jet and heavy flavor production, 
the effects of hadronization and higher order corrections 
are not yet well controlled theoretically. For detailed dis- 
cussions, see recent review paper in Ref. [4]. Therefore, 
to develop more practical ways for measuring the gluon 
polarization at RHIC experiments is one of the urgent 
theoretical tasks today. 

The quarkonium production and decays have long been 
taken as ideal processes to investigate the nature of QCD. 
Due to the approximately non-relativistic nature, the 
heavy quark and antiquark bound state is considered as 
one of the simplest system in QCD. The very clean sig- 
nals of quarkonium leptonic decays lead to the experi- 
mental observations with a high precision, and therefore 
the quarkonium may play a crucial role in investigating 
many quantities and phenomena, e.g. the parton distri- 
bution, the QGP and even new physics. 

However, it should be pointed out that although the 
heavy quarkonium physics has been investigated for 
about thirty years, theoretical description of quarkonium 
production is still premature. Conventionally, the so- 
called color-singlet model (CSM) was widely employed 
in the study of heavy quarkonium production and decays 
[5]. In CSM, it is assumed that the QQ pair produced in 
a high energy collision will bind to form a given quarko- 
nium state only if the QQ pair is created in color-singlet 
state with the same quantum numbers as the produced 
bound states; as well, in the quarkonium decays the an- 
nihilating QQ pair will be in short distance and singlet 
with the same quantum numbers as its parent bound 
states. It is assumed in CSM that the production am- 
plitudes can be factorized into short distance and long 
distance parts. The short distance sector is perturbative 
QCD applicable, while all the long distance nonpertur- 
bative effects are attributed to a single parameter, the 
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wave function. Nevertheless, color-singlet factorization 
is only an ad hoc hypothesis. There are no general argu- 
ments to guarantee such a naive assumption to be held 
up to higher order radiative corrections. Recently, from 
the Fermilab Tevatron data collected in the 1992-1993 
run, CDF group found [6-8] that the large-pT J/tpiip') 
were produced much more than the leading order(LO) 
CSM predictions. This phenomenon was referred as "■0- 
surplus" production or "-anomaly". 

In order to explain the experimental results of tp sur- 
plus production, the color-octet mechanism(COM) [9] 
was proposed based on a novel effective theory, the non- 
relativistic QCD(NRQCD) [10]. According to COM the 
produced cc pair in color-octet configuration may also 
evolve into ^ via radiation or absorption of soft glu- 
ons. For details, readers are recommended to reference 
[10]. The appearance of this new quarkonium produc- 
tion mechanism brought deep effects on the quarkonium 
physics. Having achieved the first-step success in explain- 
ing the CDF data, the color-octet mechanism encounters 
also difficulties in other processes [11]. As discussed in 
Ref. [12], the extent of importance of color-octet contri- 
butions in charmonium production remains unfixed for 
the time being. Therefore, to use quarkonium as probes 
to investigate new physics, one should really take the 
large ambiguities remaining in color octet matrix ele- 
ments and higher order corrections into consideration, 
especilly for charmonium system. 

In literature, a series of efforts have been made so 
far for measuring the polarized parton distributions 
through quarkonium production processes [13-17]. Un- 
fortunately, most of previous investigations were not di- 
rectly applied to the RHIC experiments and are spoiled 
by the uncertainties aforementioned. In this work, we 
show that the exclusive double heavy quarkonium pro- 
duction in the polarized proton-proton collision would 
be one of the ideal processes to measure the polarized 
gluon distributions and may play at least a supplemental 
role to the presently proposed program at RHIC. Due 
to being an exclusive process that can not be realized in 
COM, the uncertainties remaining in color-octet matrix 
elements are ruled out for our concern, which enables us 
to make more precise predictions. In addition, the dou- 
ble quarkonium production also has several advantages 
in reducing theoretical uncertainties from other sources. 
That is, (1) The concerned quantity, the asymmetry, is al- 
most independent of the relativistic corrections and non- 
perturbative uncertainties, which is normally not negli- 
gible in the study of charmonium system. The reason 
is that these kinds of uncertainties appearing in both 
numerator and denominator in the definition of asymme- 
try may cancel out. (2) The higher order QCD correc- 
tions can be well controlled by applying a suitable pr cut 
for the charmonium system, because the the strength of 
strong interaction here highly correlates to the transverse 



momenta of final states. (3) At RHIC energy, especially 
for p — p collision, the qq annihilation process is reason- 
ably negligible compared to the prevailing gluon-gluon 
fusion process. Therefore, the double-quarkonium stands 
as a very sensitive system for measuring gluon polariza- 
tions. 

Years ago, there were discussions of observing the 
gluon polarization by means of exclusive double J/ "0 pro- 
duction [18,19]. However, the main concerns there were 
not on RHIC physics, and analytical expressions were not 
presented for more extensive use. 
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FIG. 1. Typical Feynman diagrams for g + g ^ J/tp + J/^. 
The topological groups (A) and (B) show up in the pho- 
ton-photon to double charmonium production process. The 
group (C) is necessary in guaranteeing the gauge invariance. 

At the partonic level, the concerned process is similar 
to the photon-photon scattering in double quarkonium 
production [20]. Because of the non-abelian nature of 
QCD, here we should consider more diagrams as show 
in Figure 1. When calculating the helicity-dependent 
matrix elements, one must project onto definite helic- 
ity states for the initial giuons. In our calculation, this is 
achieved by taking the ghost free expression for the gluon 
state with helicity A = ± , 
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The double spin asymmetry A for J/ip pair production 
is defined as. 



A = 



da{p+p+ J/TpJ/Tp) — da{p+p^ — > J /ipj/ip) 



dcr{p+p+ J/ipJ/ip) + da{p+p^ 
EdAa/d^p 
Eda/d^p ' 



J/VJ/V') 
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where p+ and p- denote the helicity states of the incident 
protons being positive and negative, respectively. Here 
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and in the following, the parity conservation is taken for 
granted. In terms of the gluon densities and the partonic 
cross sections, this asymmetry reads, 

J dxidx2dAaAGixi,Q^)AGix2,Q^) 

J dxidx2daG{xi,Q'^)G{x2,Q'^) ' ^ ' 

where AG{x,Q^) = G+{x,Q^) - G_{x,Q^) and 
G{x,Q^) = G+{x,Q'^) + G-{x,Q^) are the polarized 
and unpolarized gluon distributions defined at the scale 
Q^. The unpolarized (polarized) partonic cross section a 
(Aa) is defined as, 

a = 1 ^ a (A, A'), Aa=\^XX'a (A, A'). (4) 

A,A' \,X' 

The concerned process, as schematically presented in Fig- 
ure 1, gives the polarized partonic differential cross sec- 
tion as 

dAa -16a^7r|i?(0)|'^ r„„,, ,4 

dt 81s8(m2 - t)4(TO2 - m)4 [ 

+ m^°(32110t^ + 90076iu + 32110m2) 

- 16TO^*(2025t^ + 12673t\ + 12673iw^ + 2025u^) 

+ 2t^u*(349t^ - 908t^u + lilAt^u^ - 908tM^ + 349m^) 

+ 4mi^(3903i'' + 57292*^^ + 117766*^^ + 57292^^ 

+ 3903m'') - 4m"(510f^ + 36713t''M + 135685tV 

+ 135685i\^ + 36713tu^ + 510m^) 

+ m^2(_ 1461^6 ^ 58600t^u + 364313t^M^ 

+ 594840t^u3 + 364313^ V + SSeOOiu^ - 1461u'') 

+ 4TO^t^u^(9t^ - hOhfu + AAt'u^ - 556iV 

- 556t^u^ + 44*2^5 - 505fw^ + 9u^) 

+ 2m^°(381t^ - 7111t^M - 83783^^1*2 _ i80639fV 

- 180639i^«^ - 83783tV - 7111*^*^ + 381tt^) (5) 
+ m^(-79i^ + 1272t^u + M^2Qt^u^ + 156224i^u^ 

+ 163850t^u'' + 156224^3^5 

+ 54526t\^ + 1272iu^ - 79u^) 

+ mHu{-mf- + 1471i^it + 9764iS2 _^ 12863*^^3 

+ 7196f^M^ + 12863f\^ + 9764t\^ 

+ lAlltu' - 36m^) 

- 2m''(2i'^ + m\ + 5151i^u2 + 2^'dA7t^u^ 
+ 24439*^^"^ + 24439t''u^ + 25947i3M'' 

+ SlSli^u'' + ntu^ + 2u^) , 

where i?(0) is the wave function of J/ip at the origin and 
m is the rest mass of J/ip; s, t and u are the Mandelstam 
variables for the partonic system. The sum of our results 
for different helicity combinations gives the unpolarized 
partonic differential cross section, a in (4), which agrees 
with the results in Refs. [21,22]. 



In the following, we investigate the spin asymmetry in 
the RHIC-Spin experiment at \/s = 500 GeV which is 
the planed highest RHIC energy. In our numerical cal- 
culations, the scale of the parton distribution func- 
tion and the strong coupling constant is taken to be the 
transverse momentum of J/ip for the px distributions. 
Whereas, in the calcidation of the angular distribution of 
the spin asymmetries and the integrated cross sections, 
the scale is taken to be = m? . The nonrelativistic 
relation m ~ 2 rric, with rric = 1.5 GeV, is used and 
|i?(0)|2 = 0.8 GeV^ 
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FIG. 2. Spin asymmetry of J/tp pair production in 
color-singlct model versus transverse momentum for differ- 
ent sot of GS polarized parton distributions at RHIC with 
colliding energy \/s = 500 GeV. 

We plot in Figure 2 the double spin asymmetry versus 

the transverse momentum of J/ ijj with respect to the pro- 
ton beams with different sets of the GS parameterizations 
[2] , and in Figure 3 the angular distribution of the asym- 
metry in the parton center-of-mass frame. The charmo- 
nium pair exclusive production happens only through the 
color-singlet scheme as shown in Figure 1 and so it is fea- 
sible to reconstruct the parton center-of-mass system in 
experiment. To be consistent with the use of the GS po- 
larized distributions, we use the MRST parameterization 
[23] for the unpolarized gluon distribution. From results 
shown in Figures 2 and 3, it is clear that the asymmetries 
obtained with different parameterizations diverse quite 
much. 

We notice that to observe the spin asymmetry both 
in pt and angular distributions, a high luminosity is re- 
quired. Nevertheless, having a relatively high accuracy in 
theory, even with a few experimental events our results 
may tell something about the gluon polarization inside 
the proton. From diagrams 2 and 3 it is easy to figure 
out that the charmonium pairs are produced mainly in 
low Pt and the forward direction relative to the beams, 
where the set C gives a small asymmetry compared to 
sets A and B. 
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FIG. 3. Angular differential asymmetry distribution of the 
J/il> pair in the parton center-of-mass frame at RHIC with 
colliding energy ^/s = 500 GeV. 

In order to estimate the event rate for the J/ ^ pair pro- 
ductions, we calculate the total cross section with differ- 
ent parton parameterizations. The results are collected 
in table I. As expected, the discrepancies among these 
predictions are not very large. Here, the notation cr^+^- 
mcans that the branching ratio of — > /i+zx") = 
0.0588, as the practical measuring mode to reconstruct 
the charmonium state, is included. Prom the predicted 
cross sections, wc sec that with the integrated luminos- 
ity of 800 pb^^ in the future run of RHIC, there will be 
thousands of J/tf) pair events to be detected, which can 
certainly give us some information on the gluon polariza- 
tion in the nucleon. 

TABLE I. Total cross sections for J/il) pair production at 
RHIC with ^/s = 500 GeV, evaluated with different parton 
distributions. 



Pt cut. Wc have also discussed the feasibility of observing 

the double spin asymmetry via this process. Our results 
show that for the time being the accumulated data with 
low achievable fraction of polarization in the RHIC beam 
are not yet enough to analyze this process for the purpose 
of measuring the giuon helicity distributions. However, in 
the future run with colliding energy of 500 GeV and the 
accumulated luminosity 800 pb^^, the J/tjj pair events 
can be surely detected and the gluon polarization could 
be measured. With the expected upgrade of RHIC in fu- 
ture, the J/iJj pair production may become a promising 
process with very clean signal and less theoretical uncer- 
tainties in uncovering the nucleon spin structures. 

Finally, it should be mentioned that there are two 
kinds of backgrounds which may interfere with the mea- 
surement of our proposed process. The first is the con- 
tribution of higher excited quarkonium states fccddown 
to double J/tp. Which is not negligibly small generally 
speaking, however, for pair production they are doubly 
suppressed. The second is the quakonium pairs produc- 
tion through color-octet mechanism, which is known to 
be the dominant J/ip pair production scheme at high en- 
ergy and large pT via gg gg hard interaction and with 
both final state gluons fragmenting to color-octet inter- 
mediate states and then evolving into qaurkonia nonper- 
turbativcly. Nevertheless, in both cases the J/iJj pairs 
are not exclusively produced, hence, experimentally they 
can be excluded with care. 
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To conclude, in this work we have shown that the ex- 
clusive J/tp (quarkonium) pair production at RHIC may 
stand as a novel process in measuring the gluon spin dis- 
tributions inside the polarized nucleon. The large mass 
of heavy quark guarantees that perturbative calculation 
is applicable to this process; the asymmetry, rather than 
cross sections, eliminates large amount of uncertainties 
which come from the non-perturbative hadronization and 
relativistic corrections. The higher order QCD correc- 
tions may be properly controlled by employing a suitable 
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